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Abstract. We discuss the application of the determinantal method to the 
proof of the Riemann Hypothesis. We start from the fact that, if a certain 
doubly infinite set of determinants are all positive, then the hypothesis is true. 
This approach extends the work of Csordas, Norfolk and Varga in 1986, and 
makes extensive use of the results described by Karlin. We have discovered and 
proved or conjectured relations in five areas of the problem, as summarized in 
the Introduction. Further effort could well lead to more progress. 



1. INTRODUCTION 

1.1. It is now over 150 years ago that the Riemann Hypothesis (RH) was first 
stated. Its proof or disproof is generally thought to be a very difficult problem, 
since many eminent mathematicians have so far been unable to solve it. 

In this document suggest that the determinantal approach described below 
allows the RH problem to be broken down into many perhaps simpler problems, 
which might facilitate a soloution of the whole. The determinantal approach aims 
at the proof of the inequalites of (|2.6p . The method is an extension of the work 
in 1986 of Csordas, Norfolk and Varga Q] ( see also Csordas and Varga 12]), who 
proved (|2.6[) for the case of order r = 2. We are not aware of any other significant 
progress on the determinantal approach since then. 

In many instances we use the computer to explore the problem for possible rela- 
tionships, to uncover apparent behavior of the functions involved as their variables 
change, and to carry out algebraic operations that would otherwise be extremely 
tedious. Thus we arrive at two types of result. In some cases, often stated in 
lemmas, we have rigorous statements, although often they were discovered by cal- 
culation. Alternatively sometimes we can only conjecture a general result, normally 
as an extension of patterns uncovered by calculation. In three cases we present ex- 
plicit conjectures for important proposed results, although there are others that are 
implicit. 

After a survey of previous relevant work in Sec. 2, we report new developments 
in five topics related to the determinantal method. In Sec. 3 we point out that 
the method used to prove the case of order r = 2 has a natural generalization to 
higher order provided by the sign-regularity ideas described by Karlin [3], most 
importantly the appropriate constant sign of the Wronskian of <&(u) of order r. 
Here defined in ()2.1[) . is intimately connected with the Riemann £- function. 

We refer to our recent reports [B] , [7] that use two different methods that prove the 
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validity of (|2.6|) for r = 3 and n > 2. Numerical studies make it almost certain 
that the same holds for order r — 4, except that the range of validity is then n > 3. 
These studies also make it certain that the generalized method does not work for 
r > 5. Of course it can be shown by rigorous numerical calculations that (|2.6|) does 
hold for the excluded cases for r = 3, 4, and indeed we have not found an example 
where the condiition fails in more than 100,000 tests - otherwise this report would 
have no point. 

In Sec. 4 we propose an extension of the generalized method used in Sec. 3. The 
key element in that method is a kernel K(u,v) = $(u + v). In Conjecture [T] and 
the associated discussion we propose that the cumulants {^ m (u)} of <I>(it) defined 
in (14. ip . if used in the method in place of $(«), will lead to the proof of (|2.6p 
for each r, provided that n > r)(r), where r](r) is a non-decreasing function of r. 
This conjecture is supported by numerical calculations presented in Sec. 4.4. Thus, 
even if the Conjecture were true, for all r > 2 there would be a set of values of n, 
increasing in number with r, for which the method would be unable to prove (|2.6I) . 

In Sec. 5 we discuss one such case, r = 3, n = 2, and show that D(2, 3) > will 
hold provided that q(u, v) < 0, u, v > 0, where q(u, v) is a fairly simple function of 
u, v, analogous to the Wronskian, that is defined in (|5.18|) . The condition on q(u, v) 
remains to be proved. If that can be accomplished, this example would provide 
hope that analogous methods will handle the other exceptional cases. 

In Sec. 6 we turn to an examination of the Wronskian of $(it) defined in (|6.1|) . 
The series (|2.1j) defining $(it) converges rapidly, the more so as u increases. Thus 
we can obtain a good approximation to the Wronskian for large u if we replace the 
series by its first term. That leads to the study of the determinant W r (u) given 
by (|6.3p - examples up to r — 7 are listed in (|6.4I) . From its definition it is clear 
that W r (u) is a polynomial of order r 2 , but the examples suggest that the lowest 
r(r — l)/2 coefficients are zero, and likewise for the highest r(r — l)/2 coefficients. 
The elements in (|6.3|) are polynomials {pk(y)} given by a recurrence relation [2j 
p. 184]. We have found two alternative representations for the coefficients of the 
polynomials {pk(y)}. Using the first or lower representation we have shown in 
Lemma 6.1 that the behavior observed for the lowest r(r— 1)/2 coefficients in W r (u) 
holds for all r. In Conjecture 2 we propose that the observed behavior of the highest 
r(r — l)/2 coefficients in W r (u) also applies for all r. We expect that the second or 
upper representation will be useful in the proof of this Conjecture. 

It is also clear from Sec. 2.8 and [7] that the behavior of the Wronskian (|6.ip for 
r = 2, 3 is dominated by the first term of (|2.ip . In the case of r = 3 the contribution 
from the second term in the series is significant. This contribution may also be 
expressed in terms of determinants that are polynomials [7J (3.6), p. 6]. Thus we 
see that the behavior of W r {u) and a few corresponding later determinants is crucial 
to the proof that D(n,r) > for r — 2,3,4 (except for the exceptional cases). It 
might well be that the above representations are a key factor that will lead to a 
deeper understanding of why the RH holds. 

In Sec. 7 we discuss the question of the behavior of D(n, r) as n — > oo for n 
fixed. Our original approach [4] , [5] was to use the Laplace method to approximate 
the integrals appearing in the definition of D(n,r). We found that a proof that 
D(n, r) > for large n is probably feasible provided that the truth of Conjecture 
3 can be demonstrated. This conjecture is an algebraic relation not special to the 
RH. In Sec. 7.3 we suggest an alternative approach that is also probably feasible. 
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It is based on the properties of the determinant W r (u), in particular on the re- 
quirement that the highest non-zero coefficient of W r (u) has an appropriate sign 
e r = (— l) r (' -1 )/ 2 as stated in Conjecture 2. 



2. PREVIOUS RESULTS 

2.1. The foundations of the determinantal approach may be found in several 
sources, perhaps most importantly the book of Karlin .3 . A key function in 
the notation of p] and [2], is defined by 

oo 

(2.1) $(it) = (2m 4 7r 2 e 9M - 3m 2 ne 5u ) exp (-m 2 7re 4 "). 

m— 1 

The function $(u) is even and analytic for all u with \^su\ < 7r/4 (see pQ. 
Let the Fourier cosine transform of be 

p oo 

(2.2) E(t) = / du<S>(u)cos(ut). 

Jo 

It is well known that the RH is equivalent to the statement that all the zeros of 
S(t) are real. 

In 1925 Polya therefore raised the question of what properties of the function 
<I>(it) are sufficient to lead to a proof that all the zeros of S(i) are real, but he 
presented no answers. To investigate this question we set z — —t 2 and F(z) = 2(t), 
so that equation (|2.2p leads to 



(2.3) E(t)=F(z) = J2^ zn - 



Here the series coefficients are 

(2.4) p n = 1 I du^(u)u 2n , n = 0,l,. 

1 (2n + 1) 



We define the matrix 



Pj-i, 3 > i; 



(2.5) B(i,j) = i ■>.-<{■ = 0,1,2,.... 

As Karlin [3l p. 393, p. 412] explains, the RH requirement that all the zeros of F(z) 
are real and negative is equivalent to the conditions 

(2.6) D(n,r)>0, n = 0, 1, . . . ; r = 1, 2..., 
where we denote the minors D(n, r) of order r by 

(2.7) D(n, r) = det[B(i, j + n)]i,j=i,..., r . 



The proof of the inequalities (|2.6|) is the set of simpler problems referred to above, 
and we call this approach the determinantal method. 
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2.2. In the case of r = 1 in (|2.6|) the result is trivially correct, as all coefficients 
(|2.4j) are positive, since $(it) > for all real u. 
For r = 2 we have for all n > 



(2.8) D(n,2) 



Pn Pn+1 
Pn— 1 Pn 

so that the condition (|2.6[) means that 

(2.9) 02>A,-iAh-i> n = 0,1,2,... 

The proof of (|2.9j) may be regarded as the first step in the program we are suggest- 
ing. 

In 1986 Csordas, Norfolk and Varga [2] (see also Csordas and Varga [2]) proved 
the validity of the Turan inequalities (postulated about 50 years earlier by Polya), 
which may be written as 

(2.10) 0* > Vl±±lp n _ 1 p n+1 , n =l,2,... 

n 

Since > 1 it follows that (|2.9[) is correct for n = 1, 2, . . ., and the case n = 

is obvious. Thus the proof of the first problem (|2.9j) is implicit in the results of 
Csordas, Norfolk and Varga. 



2.3. We now outline a proof of (12.91) using the basic ideas of [T] , but which is 
simpler than their method since the result is weaker, but adequate for our purposes. 
We define the kernel 

(2.11) K(u,v) = $(u + v), < u,v 

Suppose that 



(2.12) A(t) = / dv<j>(v,i)<Z>(v), 

Jo 

with the kernel 4>{v, t) being 

(2.13) 0(„, t ) = |!_, t >0. 



Note from (|23|l that 

(2.14) (3 n = X(2n + 1) n > 0. 
Let the kernel A(s,t) be defined by 

(2.15) A(s,t) = X(s + t), s,t>0. 

The argument underlying the development of [3, p. 140] shows that, for appropriate 
values of s. t, 



(2.16) J du I dv4>(u,s)$(u + v)4>{v,t) = / dv<f>(v, s + 
Thus we obtain 

(2.17) A(s,t)=l du I dv(j)(u : s)K(u,v)(j)(v,t). 
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2.4. We make considerable use of the concepts of sign-regularity and the compound 
kernel. Following Karlin [3] suppose that we have a kernel Q(x,y). Let X be a 
linearly ordered set, discrete or continuous. For a given positive integer p the open 
simplex A p (A) is 



(2.18) 



A p (X) = {x — {x\,X2, . . . , x p ) \xi < x 2 < . ■ . < x p : x, t e X} , 



and similarly for Y, y. 

As in |31 p. 12] the compound kernel Q\p\ (x,y) is defined by 



(2.19) 



Q[ P ] (x,y) 



Q(xi,vi) 



Q{X 2 ,V2) 



1 (x p ,yi) Q(x p ,y 2 ) 



Q{xi,y P ) 
Q (x2,y P ) 

Q {.Xp-, y P ) 



Suppose that the sequence e p = (— 1) 



p(p-l)/2 



p = 1, . . . , r. Then we say that 



Q(x, y) is sign-reverse regular of order r (i.e. RR r ) if e p Q[p\ (x, y) is a non-negative 
function on A p (X) x A p (Y) for each p = 1,2,..., r. Similarly, if e p = 1, p = 
1, . . . , r, then Q{x, y) is totally positive (i.e. TP r ). If the subscript r is omitted in 
the designation we mean that the sign-regularity applies to all values of r. 

From two applications of the basic composition formula (BCF) [3j (2.5), p. 17] 
to equation (I2.17[) we find that for any p > 



o Jo 



(2.20) A M (s,f) 
where u, y, s, t are defined as in (|2.18[) . 



2.5. The relevance of the above relations becomes clear when we combine the 
information in (I2.20|) . (|2.15j) and (|2.14l) . Suppose in (|2.20p we choose p = 2, and s, 
t such that 



(2.21) 



si = ii = n — 



so = <2 = n ■ 



n = 1,2. 



3 

2' "* "* " ' 2 
The elements of the determinant Q(n) = A[ 2 ] (s,t) are then 



(2.22) fi(n) = 
Thus, using ()2.14|) . we have 

(2.23) fi(n) = 



A(si+ti) A(si+< 2 ) 
A(s 2 + ti) A(s 2 + t 2 ) 



, n = l,2,. 



Pn-l Pn 



1,2,.. 



2.6. It may be shown that the kernel <f>(v,t) is TP, i.e. totally positive, which 
means that the determinant </>[ p ] (u, s) is positive for (u, s) in the appropriate simplex 
product for any p = 1,2,..., and in particular for p — 2. 

Now suppose that KnAu, v) < for all values of the arguments (i.e. K(u,v) 
is i?i? 2 ), then it follows from (|2.20j) that A(u,v) is RR 2 , so that fi(ra) < 0, n = 
1,2,.... We observe that, for n > 1, the determinant D(n, 2) = — 51(n), since one 
determinant is the other with rows reversed. 

The conclusion is that D(n, 2) > 0, n = 0, 1, . . ., as required by the RH, provided 
that K(u, v) = $(u + v) is RR 2 . 
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2.7. An efficient method of testing for sign-regularity in the present case is de- 
scribed by Karlin [3]. We have 

Lemma 2.1. Suppose that ip(x) is analytic in a neighborhood of X = Y = (0, oo), 
and that the kernel K(x,y) = i\){x + y), with x,y € (0,oo). Define w p (u) = 
det |i/> (i+J ~ 2) 0)|*\ = r If e p w p (u)>0 1 u>0, p=l,2,...,r thenK(x,y)isRR r . 

Proof. This result is a special case of [3l Theorem 2.6, p. 55]. The analyticity of 
K(x, y) ensures that the differentiability requirements of the theorem are satisfied. 
The relation 

Qi+j-2 V 



(2.24) 



det 



dx l ~ 1 dy 



— K i x ,y) 



det 



(i+J-2) 



(»..) 



*>i=i 



together with [3J (1-3), p. 48] demonstrates that the requirements on the compound 
kernel appearing in the statement of the theorem hold. □ 

Remark 2.1. In the case of r = 2 studied by [1] and [2], an alternative terminology 
is used. Instead of stating that ^(x+y) is RR2 they require that log ip(x) be strictly 
concave on (0, 00), i.e. that the second derivative of \ogi)(x) be negative. The two 
terms both lead to the condition (|2.25[) . 

To apply this lemma to the kernel K(u,v) of (|2.11[) we set ip(x) = $(2) and 
choose r — 2. The required analyticity follows from [TJ Theorem A, p. 523]. 

For p = 1 the condition of (I2.24[) requires that tp(u) > 0, u > 0, since t\ = 1 . 
This is true since $(u) > 0, u > 0. 

For p = 2 the condition becomes 

$W(u) 



(2.25) 



$(u) 
$W(u) 



$(a)( t 



$ (1) (u) 2 - $(u)$ (2) (m) > 0, it > 



since ei = — 1. 

In their discussion of the Turan inequalities, Csordas and Varga [2] study the 
case when ip(x) = &(y/x). After "a series of involved, but elementary, estimates" 
they prove that 

(2.26) u[$ (1) (u) 2 -$(u)$ (2) (u)] + $(u)$ (1) (m) >0, u>0 
so that 

(2.27) $ (1) (u) 2 - $(u)$ (2) (u) > -u _1 $(u)$ (1) (zt) > 0, u > 

Inequality (|2~2"5|) follows since $ (1) (u) < 0, it > and $W(0) = (see 0). 

Consequently we see that the work of [2] shows that K(u, v) — $(u + v) is i?i?2- 

2.8. Here we present a self-contained, simpler proof of (|2.25|) . Using the notation 
of [1 p. 184] we have 



(2.28) 

where 
(2.29) 



$00 (t 



aC?)(«) = 



n=l 



oW^tt), j = 0, 1, . 



2 p :)+ i(7rn 2 e 4 ") exp(5u - 7r?i 2 e 4u ). 



With y = 7re 4u define 

(2.30) Qj(y) = $ (j) (M)exp[-5u + 7re 4,I ]/7r, j = 0,1,2, 
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so that (|2.25|) is equivalent to 

(2.31) W(y) = n (y)Sl 2 (y)-n 1 (y) 2 < 0, y > n. 
We may write 

(2.32) n j (y)=p J+1 (y) + T J ( y ), j = 0,1,2. 
where, as in [2] 

pi(y) = -3 + 2y 

(2.33) p2(y) = -15 + 30J/-82; 2 

Mv) = -75 + 330y - 224y 2 + 32y 3 . 

In {2jj2} we define 

OO 

(2.34) T J (y) = e^n 2 p J+1 (n 2 y)exp(-n 2 y), j = 0,1,2. 

n=2 

Thus, substituting (|2.32[) into (|2.31[) . we obtain 

5 

(2.35) W(y) = W 2 (y) + Y,Ti(y) 
where 

W 2 (y) = pi(y)p 3 (y) ~Vi{vf 
Ti{y) = Pi{y)^2{y) 
/o or\ T ^(y) = P3(y)Yo(y) 

{ ] T 3 (y) = ~2p 2 (y)T 1 (y) 

Uy) = r (y)T 2 (y) 
T 5 (y) = -Ti(y) 2 

The proof of (|2.25p is based on the following three Lemmas. 
Lemma 2.2. The function W 2 (y) < -843 for y>n. 
Proof. Using the definitions (|2.33j) it may be shown that 

(2.37) W 2 (y) = 16y(-15 + I2y - Ay 2 ) = 16y[-4(y - 3/2) 2 - 6] 

Thus W 2 (y) decreases steadily for y > 1.5, so that W 2 (y) < W 2 (n) = -843.4... if 
y > 7r and the lemma is proved. □ 

Lemma 2.3. For y > tt we have \p k {y)\ < 2 {2k ~^y k , k = 1, 2, 3. 

Proof. For k = 1 the proof is obvious. For k — 2 the zero of 

(2.38) p 2 (y) + 8y 2 = 30y - 15 

is at y = 0.5 Thus the function 30y — 15 is positive for y > tt > 0.5, so that 

(2.39) p 2 (y)>~8y 2 , y > tt 
Also we find that the zeros of 

(2.40) p 2 (y)-8y 2 = -16y 2 + 30y-15 

are both complex. Thus the function — 8y 2 + 30y — 15 is negative for y > tt > 0.5, 
so that 

(2.41) p 2 (y)<8y 2 , y>n, 
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from which the proof follows. For k = 3 the proof is similar. We find that the real 
zeros of pz{y) — 32y 3 are < 1.193, while the real zeros of P3(y) + 32y 3 are < 0.275, 
which leads to the desired result. □ 

Lemma 2.4. For y > ir and j =0,1,2 we have \Tj(y)\ < 2^+^Cj(y)y^ +1 h- 3y 
where C J (y) = (l-e- 2 y2^ +4 '>/2)- 1 . 

Proof. Csordas and Varga [U Lemma 3.1, p. 185] have shown us how to prove this 
result. Using (|2.34[) and the bounds on Pj(y) of Lemma 2.3, we find that 



(2.42) 

We may write 
(2.43) 



\Tj(y)\ < 2 2j+ y +1 e y y n 2j+4 exp(-n 2 y) 



71=2 



,23+4 



cxp(-n y) = 



exp[n(ny — (2j + 4) logn/n)] 
It is not difficult to show that, if n > 2 and 2j + 4 < ir, then 

(2.44) (ny- (2j + 4) log n/n) > log K(ri) 
where K(n) = e 2a 2-^'+ 4 )/ 2 .Thus 

OO OO 

(2.45) n 2j+4 ex P (-n 2 y) < £ K{n)- n = K^C^y) 

n=2 n=2 

and the result follows. 

Combining the above information we are led to 
Lemma 2.5. 



□ 



(2.46) 



$(*)(«) $( 2 )(u) 



< 0, u > 



Proof. Above we explained how ()2.31j) is equivalent to (|2.46j) . To prove (|2.31j) 
we use Lemmas 12.21 - 12.41 to provide upper bounds to the absolute values (where 
appropriate) of all the quantities listed in (|2.36[) . which for y > it leads to 



(2.47) 



From (|2.35|) we deduce that an upper bound for W(y) is the sum of the bounds in 
(|2.47| . which is -635.80. Thus the lemma is proved. □ 

Remark 2.1. As we find by doing precise numerical computations, this bound is not 
tight. A much closer bound could be obtained by extending the method to explicitly 
include the terms corresponding to n = 2 in (|2.32|) . Our systematic approach may 
be well suited to extension to orders r = 3, 4, as the large number of terms involved 
can be handled on a computer (but still preserving rigor). 



W 2 (y) 


< -843.19 


= -843.19 


\Ti(y)\ 


< 2 14 C 2 (7r)7r 4 e- 37r 


= 132.76 


\T2(y)\ 


< 2 10 Co(7r)7r 4 e- 3 - 


= 8.30 


\Ts(y)\ 


< 2 13 Ci(7r)7r 4 e- 3 * 


= 64.88 


\Uy)\ 


< 2 18 C (7r)C 2 (7r)7r 4 e- 6 - 


= 0.17 


\T 5 (y)\ 


< 


= 
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3. DOES THE OBVIOUS EXTENSION WORK? 



3.1. The way in which we have described the work of [T], [2] in Sec. 2 suggests 
an obvious method for attempting to extend the proof of the conditions (|2.6p to 
higher values of the order r. We can adjust the procedure of Sec. 2.5, and Lemma 
12.11 applies to all values of r. We need to determine whether the kernel K(u,v) of 
(|2.1ip has sign- regularity pattern RR r for values of r > 2. 

An easy way to attack this question is via a numerical computation. With 
ip(u) = $(?/), the dervatives of i/j( u ) f° r use m (|2.24l) may be written as rapidly 
converging series obtained from (|2.ip . which can be approximated to an accuracy 
of hundreds of decimal places with little trouble. Calculations for a set of values of 
u indicate that 

(3.1) e r w r (u) > 0, u > 0, r = l,2,...,4 

but the inequality is definitely not true for order r = 5. If our code is not in error, 
then K(u,v) is not RR5, and the method of Sec. 2 cannot be used to prove any of 
the inequalities of (|2.6p for r > 5. 

We discuss a method that may to some extent overcome this difficulty in Sec. 4, 
but next we report on the case of r = 3. 

3.2. We have produced a proof [6]) that 

(3.2) D(n,3)>0, 71 = 2,3,... 

The method generally follows the pattern of [T] and [2], in that some extensive 
computations (not numerical, but computer assisted) are used to prove that the 
kernel K(u, v) is RR3. Equation (|2.20p with p = 3 then shows that Ar 3 i(s, i) is 
negative for all (s,t) provided that the components of the arguments are all > 0. 

A more systematic method of proving the same result was described in a later 
report [7]. 

To prove (|3.2p we choose 

(3.3) , 



s 2 = t 2 = n 



s 3 = t 3 = n 



n = 2,3,. 



2 2 2 

and as in Sec. 2.5 use the relations (|2.19l) . (|2.15p . (|2.14l) . which give 

Pn~l fin 

(3.4) A [3] (a,£)= p n -i p n P n+ i <0, n = 2,3,... 

Pn Pn+l Pn+2 

The restriction n > 2 arises from the need to satisfy the condition s > —1/2 from 
(l2~T3l) . 

Reversing the order of the rows of (|3.4p , thereby changing its sign, produces the 
determinant D(n,3) (remember that £3 = —1) , so that 

fin fln+l Pn+2 

(3.5) D(n,3)= /?„_! I3 n (3 n+1 

Pn-2 fin-l fin 

Just as for order r = 2 the case n = is trivial, but the case n 
in that the relation 

Pi h 83 
£1(1,3)= A, B x 3 2 >0, 
8 81 



> 0, n = 2,3,.. 



1 is exceptional 



(3.6) 
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does not follow from the general method. Its validity may be checked by inserting 
the numerical values of Pj listed by [U p. 540]. 

In Sec. 5 below we present some further observations on the case of D(l,3). 



3.3. Suppose that our calculation leading to (|3.1[) is correct for r = 4, which means 
that K(u,v) is RR4. Then we can apply reasoning similar to that in Sec. 3.2. We 
choose 



(3.7) 8j = t, 

with the result that 



(3.8) A [4] (s,i) 



n-5/2 + 2j, j =0,1, 2, 3; 



3,4. 



Pn-3 


Pn-2 


Pn-l 


Pn 


Pn-2 


Pn-l 


Pn 


Pn+1 


Pn-l 


Pn 


Pn+1 


Pn+2 


Pn 


Pn+1 


Pn+2 


Pn+3 



<o, 



3,4,.. 



In this case reversing the order of the rows does not change the sign of the 
dererminant, but, since £4 = 1, we still obtain 



(3.9) 



Pn 


Pn+1 


Pn+2 


Pn+3 


Pn-l 


Pn 


Pn+1 


Pn+2 


Pn-2 


Pn-l 


Pn 


Pn+1 


Pn-3 


Pn-2 


Pn-l 


Pn 



>o, 



3,4, 



Now there are two non-trivial cases that are not covered by the general method, 
i.e. D(l,4) and D(2,4)), but again a numerical calculation shows that these de- 
terminants are positive. We note that, in all the exceptional cases, one or more 
elements of the determinant is zero. 

4. CUMULANTS 

4.1. It appears that the answer to the question of the previous section is that the 
obvious generalization of the method of jl] and [2] probably does allow a modest 
extension to orders 3 and 4. To proceed further we present a conjecture that could 
provide a partial proof of the RH inequalities for all orders. Even if the conjecture 
is correct, there will remain to be studied an infinite set of exceptional cases such 
as occur for orders 3,4. A beginning to this study is found in Sec. [5] 
We introduce the cumulants {^ m (u)}, where 



dv&r 



(«)• 



1,2,.. 



(4.1) 
and 

(4.2) ¥„(«) = *(«). 
There is a corresponding set of kernels 

(4.3) K (u, v; m) = * m (u + v) , 0<u,v; m = 0, 1,.... 

As in (|2.16p it may be shown that, for appropriate values of m, s, t, 

y* 00 />oo />OG 

(4.4) / du dv<l>(u,s)ty m (u + v)<f>(v,t) = / dv<j>{v,s + t)9 m (v) 
Jo Jo Jo 

Integrating by parts m times leads to 



(4.5) 



dv<j>(v,s + t)y m (v) 



dv(j){v, s + 1 + m)<S>(v) = A(s + t + m) 
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where t, s > and m = 0, 1, 2, . . .. Thus we obtain 

/>oo />oo 

(4.6) A (s + to/2, t + to/2) = / / dv<p (u, s) K (u,v;m) cj) (v,t) . 

Jo Jo 

As for l|2.20[) . using the compound kernel K m (tt, u; m) corresponding to 
i4T(w, u : m), we have 

/■oo /-oo 

(4.7) A[p] (s + mw/2, t + mw/2) = du dv_ (f>^ (u, s) (u, u; m) (f>\^ (v, t) , 

•/•• Jo 

where u, v, s, t are defined as in (|2.18p ). and w is a vector with all components 
equal to unity. 

4.2. Now we adapt the procedure of (13.31) . (|3.7p to relate the compound kernel 
appearing in the LHS of (|4.7[) to some of the determinants needed in (|2.6p . In (14.7[) 

we choose p = r > 1 , and s, t such that 

(4.8) Sj =tj=n + 2j, j = l,...,r. 

The elements of the determinant Ar r i (s + mw/2t + mw/2) are then 

(4.9) A(2^ + 2i + 2j + TO), i,j' = l,...,r 

If the argument of A appearing in (|4.9j) is odd and positive for all entries, then we 
can use (|2.15j) to express all elements of Ar r i (s + mw/2t + mw/2) in terms of the 
coefficients (3 n given by ()2.14j) . For use in (|4.7j) the inequality (|2.13|) requires that 
fi + 2 > 0. It follows that the determinant with elements given by (j4.9[) is the same 
as D(k,r) with the order of the rows reversed, and an appropriate choice of k. As 
before, the reversal of the rows changes the determinant by a factor of e r . 
It follows that all values of k are possible so long as k > k^. It is found that 

(4.10) k L =n L +r-l, 
where til — to/2, to even: (to + 1)/2, m odd. 

4.3. Let us assume the validity of the following conjecture 

CONJECTURE 1. For any given order r > 1, there is a lowest integer m(r) such 
that K(u,v;m(r)) is RR r . 

Remark 4.1. Karlin 3, (1-17), p. 102; Remark 5.6, p. 128] shows that, because of 
the relation (|4.ip . if K(u,v;mi) is RR r , then so is K(u,v;m2) for all mi > m\. 
Thus the conjecture implies that, if to > m{r) 1 then K(u,v;m) is RR r . It also 
follows that, if r-i > r*i, then m(r2) > m{r\). 

The Conjecture means that the compound kernel K< r -i (u,v; m(r)) has sign e r , so 
that, as before, the same applies to A[ r ] (s + mw/2, t + mw/2) from (|4.7I) . since the 
other factors in the integrand are non- negative. 

The conclusion is that, given the Conjecture, we can prove the positivity of 
D(n,r) needed in (|2.6|) for all values of n,r such that n > fci(TO.(r)). 
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4.4. To examine in a non-rigorous numerical way whether the conjecture might 
be valid we have applied the method of Sec. 3.1 to the kernels K(x,y;m), m — 
0,1,2,...; r = 2,3,..., with each sequence going beyond 16. The results are 
consistent with the Conjecture, some values of m(r), r)(r)) being listed in the Table, 
where rj(r) = ki,(m(r)) denotes the number of exceptional cases (not including 
n = 0) for the indicated value of r. 



r 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 
m(r) 1 1 1 2 4 6 7 7 7 9 11 13 15 15 15 16 
r)(r) 1 2 3 5 6 7 8 10 12 14 15 16 18 20 22 24 25 26 27 

If the conjecture is correct, then the proof of the RH reduces to showing that 

(4.11) D(n,r)>0, n = 0, . . . , r)(r); r = 1, 2..., 

4.5. In some applications it might be helpful to have an explicit form for the cumu- 
lant $ m (u) rather than the iterative relations (|4.1[> . Karlin [3] discussed cumulants 
(although apparently with a different aim in mind), and he showed that 

1 C°° 

(4.12) * m ( tt ) = — — / dt$(t)(t - u) m -\ m > 1. 

T(m) J u 

This formula may be verified by differentiating with respect to u, which shows that 

(4.13) = -* TO _i(t»), m >i, 
consistent with (14.111 and (14.21). 



5. THE EXCEPTIONAL CASE £>(1,3) 

5.1. In Scc.[3]we explained that the relations (|3.1[i e r w r (u) > 0, u > 0, r = 1, 2, 3 

led to the result Q3.2p that 3) > 0, n = 2, 3, In this section we describe an 

analogous technique, which shows that the relation D(l, 3) > holds provided that 
q(u,v) < for all u, v > 0, where q(u,v) is defined in (|5.18|) . This is the simplest 
example of an exceptional case. 

We modify the development of Sec. 2.2 by replacing $(it) by f(u) = ^^'(u), and 
defining the kernel F(u,v) = f(u + v). Note that [TJ p. 523] states that /(0) = 
because $(u) is an even function of u. Thus, with s, f > 0, (|2.17l) becomes 



(5.1) du dv4>(u,s)f(u + v)(f>(v,t) = / s + t)/(u) = w(s + i), 
Jo Jo Jo 

where 

/•OO 

(5.2) il(s,t)=u(s + t)= dv(/)(v,s + t)f(v). 

Jo 

As in (|2.20p we apply the BCF with p = 3, and find that 

/•oo />oo 

(5-3) fi [3] (s,t)= du dv(p [3] (u,s)F [3] (u,v)<f> [3] (v,t), 

Jo Jo 

where u, y, s, t are defined as before, including inequalities for u, y as in (|2.18l) . 
We replace the RHS of ((531) using [3j (2.6) p. 17], so that ((531) becomes 



(5.4) %](&*)=/ du / du0(u,5)F [3] (u,u)e(u,t) 
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where 

(5.5) Q(u, s) = 4>{ui, s\)<p(u2, S2)4>{m, S3) 
and all the integrals now run over (0,oo). 

5.2. Now in (|5.4[) we choose s, t with z > such that 

(5.6) s j =t j =z + 2j-2, j = l,...,3. 
The elements of the determinant firai (s, i) are then 

(5.7) A(i,j) = w(22 + 2i + 2j-4), *,j = l,...,3. 

The next step is to take the limit as z — > of both sides of (|5.3p . On the LHS 
consider first the element 

(5.8) A(l,l)=w(2z)= H dv^—Q^iv) 



T(2z) 

Integrating by parts (differentiating ^>^(v)) leads to 

(5.9) A (l,l) = — — - / dvv 2 *& 2 Hv) 

T(2z + 1) J 

so that 

(5.10) lim A(l, 1) = 0, 



since 
(5.11) 



/ dw$ (2) (v) = (00) - (0) = 
Jo 



on account of the properties of $W(it). 
For the remaining cases we have 

(5.12) A(iJ) = l dv T(2z + 2 . + 2j _ A) ^(v), ij = 1,2,3; i +j >2. 
Integrating by parts (integrating $W(i>)) and setting z — leads to 

(5.13) A(i,j) = TVo . i0 . - / dwP+*-**(v) = A(2z + 2j - 5) = -ft +J -_ 3 



r(2t + 2j - 3) y 

Thus for z — the elements of Orai (s,t) are the negative of the elements of the 
determinant D(l, 3) (see ([376]) ) with the order of the rows reversed. 

5.3. We now turn to the RHS of pT3]) . The function 1/T(z) = O(z) as 2 -> , 
so that in that limit the contribution to (|5.4p from any region of integration with 
iti > uq or wi > uo for any positive uq will be zero. Since 

/>oo 

(5.14) r(z)= / diiu'^e^, 

Jo 

we end up with 

Lemma 5.1. Apart from an inessential positive constant factor, 
(5.15) 

/■OO /" OO />OG />00 

D(l,3)=/ eht 2 / du 3 / ^2 / d« 3 ^(u 2 ,2)0(u 3 ,4)F [3] (u,w)0(«a,2)^(w 8 ,4) 
Jo Jo Jo Jo 

w/iere u = (0, u 2 ,u 3 ),v = (0,v 2 ,v 3 ) 
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Proof. In case the above remarks are not considered to be an adequate proof of this 
lemma, we outline an alternative method. Expand each of the two determinants in 
(|5.15[) into their four components (remember that the (1,1) element of Fw(u, v) = 
f(u\ + vi) = ^^'(O) = 0), and perform the integrals. It will be found that the four 
terms on each side of the equation are equal term by term. For example, consider 
the component from Fig\(u, v) of the form f (u 2 ) f (v 2 ) f (u 3 + v 3 ), which leads to 



(5.16) 



du 2 du 3 dv 2 dv 3 cj)(u 2 ,2)(j)(u 3 , 4)0(v 2 , 2)(j)(v 3 , 4) / (u 2 ) f '(v 2 ) f %u 3 + v 3 ) 



C 



= C 



du 2 u 2 f(u 2 
duu$^(u) 



dv 2 v 2 f(v 2 ) 



du 3 dv 3 u 3 3 v 3 3 f(u 3 + v 3 ) 



du / dw-u 3 w 3 $ (1) (M + v) 



The last line is easily shown to be a multiple of (3q(3 3 , where (3 n is given by (|2.4|) . 
When the appropriate constant factors are included, a similar procedure for the 
other three terms leads to Lemma 1531 after taking account of (|2.7[) . □ 

5.4. To further examine the case of D(l, 3) define Q{u 2 , u 3 , v 2 , v 3 ) as 

f(v 2 ) f(v 3 ) 

(5.17) Q(u 2 ,u 3 ,v 2 ,v 3 ) = F [3] (u,v) = f(u 2 ) f{u 2 + v 2 ) f{u 2 +v 3 ) 

f( u 3) f{u 3 + v 2 ) f(u 3 + v 3 ) 

where u = (0, u 2 , u 3 ), v = (0, v 2 , v 3 ), and /(0) = as mentioned above. In analogy 
with the function w p (u) defined in Lemma 2.1 we introduce q(u,v) as 

f(v) / ((1) («) 

(5.18) q(u,v)= f(u) f(u + v) f^(u + v) 

/W(u) fM(u + v) fW(u + v), 

which relates to the behavior of Q{u 2 , u 3 ,v 2 , v 3 ) for u 2 ps u 3l v 2 ~ v 3 . 
In analogy with Lemma 12.11 we have 

Lemma 5.2. Suppose that q(u,v) < for all u, v > 0, and also that f(u) < for 
all u > 0. Then Q{u 2 ,u 3 ,v 2 ,v 3 ) > for all u 2 > u 3 and v 2 > v 3 . 



Proof. Make the definitions 
(5.19) 



h(u, v) = 



f(u + v) 
/(«) 



(5.20) 



(5.21) 



R(u, v) 



1 



f(v) 2 
G(u 2 ,u 3 ,v) = 



dh(u,v) (1) ' 

f{v) — 5 / l >{v)h{u,v) 

ov 



f(v) 



[h{u 2 ,v) - h(u 3 ,v)] 



The following relations are easily checked. 



(5.22) 



d_ 

dv 



1 



■h(u, v) 



f(v} 



:R(u,v) 



(5.23) 
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(5.24) 



dR(u,v) 
du 



/(«) 



7 q{u,v) 



(5.25) Q(u 2 , u 3 , v 2 ,v 3 ) = f(u 2 )f(u 3 )f(v 2 )f(v 3 ) [G{u 2 ,u 3 , v 2 ) - G(u 2 , u 3 , v 3 )} 
The assumption of Lemma 15.21 together with (15.24[) shows that 

dR(u, v) 



(5.26) 



du 



> 



for all u, v > 0. Thus (|5T2"5|) implies that 

8G{u 2 ,u 3 ,v) 



(5.27) 



Ov 



> 



for all v > if u 2 > u 3 > 0. This in turn means with (|5.25j) that 

(5.28) Q{u 2 ,u 3 ,v 2 ,v 3 ) > 

if m 2 > u 3 and v 2 > v 3 as required in the statement of the lemma. □ 
Lemma 5.3. With the assumptions of Lemma \5.S\ then D(l,3) > 0. 
Proof. From (|5.17p it is clear that 

(5.29) Q{u 3 ,u 2l v 2 ,v 3 ) = Q(u 2 ,u 3l v 3 ,v 2 ) = -Q{u 2l u 3l v 2l v 3 ). 

Thus, using ()5.17j) for Fra (u, w), with C a positive constant, we can rewrite Lemma 
3.1 as 



o 







(5.30) D(1,3) = C du 2 du 3 dv 2 / dv 3 4>[ 2 ](u,s)Q(u 3 ,u 2 ,v 2 ,v 3 )(j)[ 2] (v 1 t) , 



where u — (u 2 ,u 3 ), v — (v 2 ,v 3 ) and s = t= (2,4). For u 2 > U3 and v 2 > v 3 we 
have 



(5.31) 



<P{2] {u,s)<p [2] (v,t) = 



U 2 U 3 V 2 V 3 
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(«3 - u l)( v l - v l) > °. 



□ 



so that (|5.30|) shows that 73(1,3) > as required. 

Remark 5.1. The assumption in Lemma l5.2l that f(u) = ^^(u) < for u > may 
be disregarded. It was proved to be true by Wintner [5]. 

6. COEFFICIENT REPRESENTATIONS 

6.1. We described in Sec. 2.8 how the derivatives of $(u) can be written in terms 
of a set of polynomials Pk{y), k — 1,2, ... , introduced by Csordas and Varga [2j 
(3.9), p. 185] . These CV polynomials play an important role in the study of the 
determinantal method. For instance consider the Wronskian of general order r 

$(u) ... & r -^(u) 

$(!)(«) $(2)( M ) ... $W( U ) 



(6.1) 



w r (u) — 



Using the first term in the expansion of <&( fe ) we see that, for large values of u, we 
may approximate w r (u) by 

(6.2) w r (u)w — ^ ^ VF r (y), 



1(5 



where 



(6.3) W r (y) 



pi{y) &s(y) ■■■ p r (y) 
P2(y) P3(y) ■•■ Pr+i{y) 

Pr(y) Pr+l(y) ■■■ P2r-l(y) 



with y = 7re 4n . We have seen in Sec. 2.8 and [7] how the properties of W r (y) 
were crucial in proving the positivity of D(n,r) for r — 2,3. Later in Sec. 7.3 we 
shall argue that the properties of W r (y) could also form the basis of a proof that 
positivity applies for any fixed r so long as n is large enough. 

The polynomials Pj(y) of degree j are defined by recurrence relations, starting 
with pi(y) = 2y — 3, that are given in Csordas and Varga [3J p. 184]. From the 
definition (|6.3[1 it is clear that W r (y) is a polynomial of degree at most r 2 . 

Using e r = (—l)^ r \ where /i(r) = r(r — l)/2, we find that 

t 2 W 2 {y) = 240y - 192y 2 + 64y 3 

e 3 W 3 ly) = -860160y 3 + 737280y 4 - 294912y 5 + 65536y 6 
e 4 W 4 (2/) = 190253629440?/ 6 - 169114337280y 7 + 72477573120y 8 

-19327352832y 9 + 3221225472y 10 
e 5 W 5 (y) = -0.329167393077068 10 19 y 10 + 0.299243084615516 10 19 ?/ 11 
-0.132996926495785 10 19 y 12 + 0.379991218559386 10 18 y 13 
-0.759982437118771 10 17 y 14 + 0.101330991615836 10 17 y 15 

(6.4) e 6 W 6 (y) = 0.538444964246560 10 28 y 15 - 0.497026120842978 10 28 y 16 

+0.225920964019536 10 28 y 17 - 0.669395448946 772 10 27 y 18 
+0.14344188191 7165 10 27 y 19 - 0.229507011067465 10 26 y 20 
+0.255007790074961 10 25 y 21 
e 7 W 7 (y) = -0.975629606681896 10 39 y 21 + 0.910587632903103 10 39 y 22 
-0.420271215186047 10 39 y 23 + 0.127354913692742 10 39 y 24 
-0.283010919317204 10 38 y 25 + 0.485161575972349 10 37 y 26 
-0.646882101296465 10 36 y 27 + 0.616078191 710919 10 35 y 28 

Note that all the coefficients in (|6.3[) are integers, but we have displayed only the 
15 leading digits. 

The examples in (|6.4j) suggest three general properties of the polynomials W r (y). 
We have proved one of these properties, which is stated in Lemma 16.11 Some 
progress has been made on proving the other two properties, which are described 
in Conjecture [2] These properties all relate to the coefficients {^{j, r)} of the 
polynomials {W r (y)} defined by 

r 2 

(6.5) e r W r {y) = Y,<3,r)y j . 

3=0 

We have 

Lemma 6.1. For j = 0, . . . , fi(r) — 1 the coefficients r) — 0, r > 2. 

Proof. The proof, given in Sec. 6.3, depends on Lemma l6.2l . □ 



CONJECTURE 2. For a given r > 1 the highest non-zero coefficient is 7(/i(r + 
1), r), and j(n(r + 1), r) > 0. 
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6.2. The above lemma depends on the existence of expressions for the coefficients 
of the CV polynomials [2j p. 184], which are defined by 

k 

(6.6) p k (y) = J2 d ^ k )y j - 

3=0 

There are two alternative representations of the coefficients {d(j, k)}, the first of 
which is described in this section. We call it the 'lower representation' because 
it is useful in proving Lemma 16.11 We expect that the second, called the 'upper 
representation', will be useful in proving Conjecture [2] It is described in Sec. 6.4. 

Lemma 6.2. For k > 1 the coefficients d(j, k) are given by 

3 

(6.7) d(j,k) = J2<J,v)(^ + 5) k -\ j<k, 

T/=0 

where 

(6.8) cM= tspmii, „<,. 



Also, with d(j, k) defined by 

(6.9) d(J,k)=0, j>k. 
Proof. The relation given by [H p. 184] is 

(6.10) Pk+1 (y) = 4yp£\y) + (5 - 4y) Pk (y), fc-1, 2,. 



Substituting (|6.6j) for k > 1 leads to 

d(0,fc + l) = 5d(0,fc) 

(6.11) d(j, fc + 1) = -4d(j - 1, k) + [4j + 5]d(j, fc), j = 1, . . . , fc 

+ 1, fc + 1) = -4d(k,k), 

where the recurrence starts with 

(6.12) d(0,l) = -3; d(l,l) = 2. 

It is obvious that there is a unique solution for d(j, fc), j < k obtained by iterating 
(16.111) on fc and starting from the conditions (|6.12l) . Therefore the lemma will be 
proved if we can show that the definition (|6.7[) satisfies (|6.11[) and (|6.12l) . 

For the first line of (|6.1ip . the formula (16.71) shows that 

(6.13) d(0,k) = -SxS^ 1 
so that 

(6.14) d(0, k + 1) = 5d(0, k) 
as required. 

With a given value of j the difference J between the two sides of the second line 
of equation (|6. 1 1|) is 

3 3-1 

(6.15) J = c(j, 7/)(4r/ + 5) fc_1 [4?7 + 5 - 4j - 5] + 4^ c(j - 1, r/)(4r/ + 5) fc_1 . 

?7=0 77=0 

Since (j — rj)c(j, rj) = c(j — 1, 77), 77 = 0, ... ,j — 1 it follows that J = as required. 
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To examine the third line we note that the definition (|6.7|) leads to 

( , m] d(k, k) = E; =q (-l) T7+1 (2r ? + 3)(4tj + $)^[ V \[k - r;)!]" 1 

' J =2x4 fe - 1 EU(- 1 )' ,+1 ^+^)] W-(k-v)r\ 

where P(rj) is a polynomial in ?y of degree k — 1. Now we define 

(6-17) /W- fc! E^^)T = (l + ^- 

Therefore, with j < k, 

(6.18) (.|)'/WI— ""El^i-HM-D 

since /^(-l) = 0, j < fc. It follows that 
Using this result in (|6 . 1 6|) shows that 

(6.20) d(k, k) = (-l) k+1 2x4 k -\k = 1,2,... 

since the contribution of the term containing P(rf) will be zero. Together the above 
results prove (|6.7|) of the lemma. Equation (|6.9j) follows in the same way using 
([615)1 . □ 

Some values of c(J, 77) are shown in (|6.21[) 



3 


r\ = 


77 = 1 


77 = 2 


^ = ; 





-3 








1 


-3 


5 






2 


-3/2 


5 


-7/2 




3 


-3/6 


5/2 


-7/2 


9/6 


4 


-3/24 


5/6 


-7/4 


9/6 



(6.21) 



-11/24 

6.3. We now use the representation of the coefficients d(j,r) demonstrated in 
Lemma 16.21 to prove Lemma 16.11 Define a set of column vectors ^(77) of dimension 
r by v(i]) = 1 (?7), Vi(ri), • ■ • ,v r (r])] T where 

(6.22) v a {r)) = (4?7 + 5) Q ~ 1 , a = l,...,r. 

Also let the column vector V((3, y) = \pp(y),Pp+i(y), ■ ■ .p/3 +r -i(y)] T , /3 = l,...,r, 
so that 

(6.23) W r {y)=det[V{P,y)]p =1 ,..., r 
Then we have, using (16.61) and (16. 7|) . 

H-/3-1 

(6.24) V{P,V)= E H{j^)y\ = 1, 2, . . . , r, 
where 

(6.25) ^(j,/3) = ^c(j,r / )(4r / + 5)^ 1 7;('7), £ = 1, 2, . . . , r. 

J7=0 
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To prove Lemma l6.1l wc note that we have shown in Lemma 16 . 21 that each column 
YX/3,y) of the determinant W r is a polynomial in y of degree r + (3 — 1. Each 
coefficient in these polynomials is a linear combination of column vectors v_{rf). 
Note that rj<j for the v(rf) vectors appearing in the coefficient of y° . 

The determinant W r {y) may be expressed as a sum of the determinants of all 
combinations of r column vectors vjrf) obtained by choosing one v(j]) from each 
column V_(/3, y), /3 = 1, . . . , r. In this process, we consider two v(r]) with the same 
value for 77 to be different if they appear in terms with different j. 

Every example of the resulting determinants will include as a factor det {v(jl(f3))] 
where ry(/3) is the value of rj chosen for the term in column /3. If for any two different 
values j3 — /?i,/32 it happens that ry(/3i) = rjifa), then the determinant will be zero 
since it contains two identical columns v{r)(/3i)) and vJeiifc))- 

If no pair of r/(/3i), r/(/?2) are equal then 



(6.26) ^ *?(/?)> ^=(r(r-l))/2 = /i(r) 

0=1 

and conversely, if 



i=0 



(6.27) ^77(^)< M (r) 

0=1 

then must be such an equal pair. If j([3) is the value j associated with rj(f3) it 
follows, since i(/?)>7?(/?), that there is an equal pair if 

r 

(6.28) ^j(/?)<m, 

0=1 



thus proving the claim of Lemma 16.11 

Remark 6.3. Using similar arguments it may be proved that the coefficient 7(/x(r), r) 
defined in (|6.5|) is given by 



3.29) 7 (m,r) = jf] j) j [det |(5 + 4») 



I ijj— 0, . . . ,r- 

This result may be extended to a few higher coefficients without much difficulty. 

6.4. We now turn to the second of the alternative representations for the coeffi- 
cients of the CV polynomials, the 'upper representation'. For k = 5, 9 we define an 
integer function s(i,j), i = 1, 2, ... by the initial conditions 

s(i,j) =0, i<0, j = l,2,... 
s(i,j) =0, j < 0, i = 1,2,... 

(6.30) s(i,0) =1, i = 1,2,... 

i = i,2,... 

«(<,») =2x(-4)- 1 , t = l,2,... 

and then recursively by 

(6.31) 8(i + l,i)=«(i,i) + a(t + l,j-l)(fc + 4t), i = 1,2,..., j = l,2,... 
We obtain some properties of {s(i, j)} in 
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Lemma 6.4. Fori,j > 1 letv_ = (i/i,z/2, ■ • ■ where {v a } are integers such that 
0<v a <j and Y,a=i v a=j- For given {i, j} let v(m), m = 1, 2, . . . , N(i,j) 
be the list of all possible partitions with i > 1, j > 0, thus defining N(i,j). Then 

(6.32) nm= Wj^ 

and 

(6.33) fl (*,j)= e n( fc +( 4 ( a - i ))^ (m) • 

m=l \a— 1 / 

Proof. The proof is by induction. First we note that, for j = 1, there are i partitions 
with v a (m) = 1, m = a; = 0, m ^ a for m, a = 1, 2, . . . , i. Thus A(i, 1) = i, in 
agreement with (|6.32p . and ()6.33j) becomes 

i 

(6.34) 8 (»,l)=£(fc + 4(m-l)), 
which leads to 

(6.35) s(« + 1,1)- s(i, l) = k + U, i = 1, 2, . . . 

This is consistent with (|6.30[) for j = 1, bearing in mind that s(«,0) = 1. 

Similarly, when i = 1, there is just one partition, which ^i(l) = j. In this 
case (|6.32p gives N(l,j) = 1, and (16.33)) leads to s(l,j) — both of which are 
consistent with (|6.30l) . 

Next we investigate the sequence {s(i, 2), i = 2,3,...}. Consider first i = 2. 
From (|6.3ip we see that s(2, 2) should be related to s(l, 2) and s(2, 1). We know that 
N(l, 2) = 1 and that the single partition for s(l, 2) is = (2), with s(l, 2) = k 2 . 
Also 7V(2, 1) = 2 and that the two partitions for s(2, 1) are v(l) = (1,0), z^(2) = 
(0,1), with s(2,l) = k+ [fc + 4]. 

According to the definition in the statement of the lemma, the partitions corre- 
sponding to s(2, 2) must have the form vjm) — (i/±(jn), U2(m)) where 

(6.36) 0<j/ q (to)<2, a =1,2; : m = 1, 2, . . . , N(2, 2) 
and 

(6.37) !/i(to) + i/ 2 (to) =2, m= 1, 2, . . . , A(2, 2). 
This means that A^(2, 2) = 3 and that the 3 partitions must be 

(6.38) 2/(1) - (2, 0), i/(2) - (1, 1), i/(3) = (0, 2) 
with 

(6.39) s(2, 2) = [/c 2 ] + [k(k + 4)] + [(fc + 4) 2 = [fc 2 ] + [k + (k + 4)] (k + 4) 

We see that (|6.39|) corresponds precisely to (|6.31j) with i = 1, j = 2, which reads 

(6.40) s(2, 2) = s(l, 2) + s(2, l)(fc + 4) 

In general s(i + 1, j) is the sum of all monomials of degree j constructed from 
products of numbers in the sequence fc, k + 4, k + 8, . . . , k + Ai. The equation 
(|6.31j) declares that this sum is obtained by taking the sum all the monomials in 
s(i,j), which are also of degree j but include only those formed from the sequence 
k, k + 4, k + 8, . . . , k + 4(i — 1). To that must be added the sum of all monomials of 
degree j — 1 in s(i + 1, j — 1) constructed from products of numbers in the sequence 



21 



k,k + A,k + 8, . . . ,k + Ai, which are then mulitiplied by k + Ai. All the monomials 
in s(i + 1, j) will appear from one or other of these two terms, and there will be no 
repeats. 

In a similar fashion we obtain 

(6.41) N(i + l,j) = N(i,j) + N(i + l,j-l) 

which is satisfied by (|6.32p . so that the Lemma is proved. □ 

We are now in a position to state and prove the CV polynomial coefficients 
'upper representation'. Let Ss(i,j),ag(i,j) denote the function s(i,j) of (|6.30l) . 
(|6.3ip . with the subscript being the value of k. 

Lemma 6.5. The coefficients {d(j,k)} defined in 16. 6]) are given by 

( , d9 , d(i,i + j) =-3x(-4)V*+U-l) 

1 ' +2 x (-i)«-V8 9 (i, j), i = 0,1,...; j = 0,1,... 

Proof. For the special cases where {i — 0, j = 1, 2, . . .} or {i = j, j = 1,2,.. .}, 
(|6.42jl may be checked explicitly. For the remaining cases we substitute the expres- 
sion (|6.42p into the recurrence relation (|6.11l) (middle line). The difference between 
left and right hand sides is 

d{i, i+j + l) + Ad(i - 1, i + j) - [4i + 5]d(i, i+j) = 

(6.43) -3 x (-4) 1 (s 5 (i + l,j) - s 5 (ij) - [Ai + 5}s 5 (i + 1J- 1)) 
+2 x {-A) 1 - 1 {s 9 {i,j + 1) - s 9 (i - l,j + 1) - [Ai + 5}s 9 (i,j)) , 

which is zero on account of (|6.31|) . thus concluding the proof. □ 

6.5. Another interesting and useful problem of a similar nature to that considered 
in Sec. 6.1 is discussed in this section. Equation (|2.4p reads 

1 f°° 

(6.44) (3 n = — — / du<A>{u)u 2n , n = 0,l,...,. 

r(2n + l) Jo 

Suppose we define 

poo 

(6.45) b n = / du<S>(u)u 2n , n = 0, 1, . . . , . 

Jo 

Then from (|2.4p and (|2.7p . with n > r, we may write 

(6.46) D(n, r) = det [ . 

. r ( 2 U -* + n) + 1). 

It is therefore of interest to examine the properties of 

(6.47) A(n,r)=det T(2n+1 



r(2(j - * + n) + 1) 



Appendix A of [1] and some results in [31 p. 106] show that A(n, r) > 0, n > r. 
Moreover, define 

(6.48) A(n,r) = (j[ 9(h A(n, r), 
where 

(6.49) 0(i,n) = (2n + 2i)(2n + 2i- 1). 
Then we have 
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Lemma 6.6. The function A(n, r) is a polynomial in y = l/(2n) of degree r(r — 1) 
for n > r. With 

r(r— 1) 

(6.50) A(n,r) = ^ 
then 

(6.51) (5(«) = 0, £ = 0,...,r(r-l)/2-l. 

Proof. The proof follows a pattern analogous to that used in Lemma l6.il From the 
above definitions it follows that 

(6.52) A(n,r)ij = g(j,n - i + l)A(n,r) iJ+1 , i = l,...,r; j = 1, . . . ,r - 1, 

with A(n, r)i. r = 1, i = 1, . . . , r. 

The elements A(n, r)ij are polynomials in y of order 2(r~j), j — 1, . . . , r. Define 
a set of column vectors w(7y) of dimension r by w(?7) = [ui(r]),U2(ri), . . . ,u r (i])] T 
where 

(6.53) ^(r?) = (i) r \ i = l,...,r. 

Also let the column vector U(j, y) = [A(n, r)ij-, A(n, r)2,j, . . . , A(n, r)r,j] T , J = 

1. . . . , r, , so that 

(6.54) &(n,r)=det\U(j,y)] j ^ tT 
It may be seen that 

2O—3) 

(6.55) C£(j,y)= 51 c o,ky k u(k), j = l,...,r 

k=0 

for some set of coefficients c^fe. 

From this point the argument of Lemma 16.21 applies word for word. □ 

7. D(n, r) FOR LARGE n 

7.1. We now turn to the question of the behavior for D(n,r) for large n with r 
fixed. In [4] we proposed a method, sketched here, for treating this problem. In 
Sec. 7.3 we discuss an alternative method. 

In the current notation it was conjectured that a scaled, normalized version of 
D(n, r) is approximately 



(7.1) D(n ir )sN~y ](r) Y, P ^ 



j{r) 1- -\j(r)-i 



yf3(n) 



00, 



i=0 

where j(r) = r(r — l)/2, y = (2n) _1 and P(i) is positive for i — j(r). As n — > 00 
the quantity 1/ (y/3)— >0 slowly, so that the term i — j(r) eventually dominates. Our 
calculations show that other terms can be significant for n in the thousands and 
beyond, the more so for higher r. 

The key reason for the behavior in (|7.1I) is an algebraic result, not special to the 
case of the RH, which means that certain terms in the expansion [H (2.26)] vanish 
identically. The relative sizes of the terms in (|7.ip are related to the connection of 
the formula to the RH. In (|7.1[) we have used the Laplace method (see Olver [5]) 
to approximate the integrals in (|2.4|) . 
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In a later report [5], for the case r — 2, we described an improved version of the 
method of [3], which may be generalized to r > 2 (see Sec. 7.3). We now present 
the main steps in the new method for r = 2. Expanding (|6.46[) gives, for n > 3 say, 

/ 7 9 n nf 9 n _ ^» b n -ib n +i 

{ ' > {n ' j_ r(2n + l) 2 r(2n-l)r(2n + 3) 

With g(i, n) as defined in (|6.49l) we normalize D(n, 2) to obtain, with r = 2, 

D(n, 2) = y 2 (n^ 1 «)•-*) T(2n + l) 2 #(n, r) 
( 7 - 3 ) = y 2 (b 2 n (2n + l)(2n + 2) - 6 n _i6 n+1 ) 

= B(l)(l + y)(l + 2y)-B(2)(l-»), 

where we have written = 6 2 ; B(2) = b n -ib n+x . We call S(l), J B(2) the 

components of the determinant det [bj-i +n ] i . =1 r . 

From (|6.45[) . after symmetrizing the double integral (without affecting its value), 
we may write 

poo poo 

(7.4) B(k)= / / dtxdt^itiWt^-H^RihMM), k = l,2, 

Jo Jo 

where 

(7.5) R(l,h,t 2 ) =t\th R(2,t 1 ,t 2 ) = (t$ + 4)/2. 

Now we set ti = r(l + Xi), i — 1,2, where we will expand the integrand about 
ti = t to be chosen later. We obtain 

ti, i 2 ) = r 4 [l + 2(xi + x 2 ) + ( Xl 2 + x 2 2 ) + 4xix 2 
/ 7fi x +2(a;i 2 .x 2 + a: 2 2 .Ti) + xi 2 x 2 2 } 

[ ' R(2,h,t 2 ) -(r 4 /2)[2 + 4(x 1 +x 2 ) + 6(x 1 2 +x 2 2 ) 

+4(a; 1 3 +x 2 3 ) + (a ;i 4 + a ;2 4 )]. 

The expressions such as 1, {x\ +x 2 )/2, {x\ 2 + x 2 2 )/2, (xix 2 ) 7 (xi 3 + x 2 3 )/2, etc. in 
(I7.6|) are symmetrized monomials. We can classify them by their type (m,j), where 
m represents the degree and j = 1, . . . , np(m) is the label of a monomial of degree 
m. For example, for m = 0, np(0) = 1 and the single monomial has structure 
1. For m = 1, np(l) = 1 and the single monomial has structure x. For m = 2, 
np(2) — 2, and the two types are x 2 and xcc. For m = 3 we have 2 types cc 3 , x 2 x, 
while for m = 4 there are 3 types . We choose the normalization of 

each symmetrized monomial such that, if x% = x 2 = x, then the expression equals 
x m . We call these monomials /i(m,j), j = l,np(m). 
With these definitions we write 

4 np(m) 

(7.7) R(k,h,t 2 ) =r 4 T(m,j,k)»(m,j), k = 1,2, 

m— j — 1 

which defines the set of coefficients {T(m, j, &)} In view of f|7.4p and f|7.5p we may 
write 

2 4 np(rn) 2 

(7.8) S(n,2)=^^ £ £ e(2,fc)z(i, fc)T(m,j,% i J(m,j) 

i=0 ro=0 j'=l fc=l 

with 

(7.9) I{m,j)=T A / dtid* a *(ti)*(t 2 )t? , - a ^ B - 2 A 1 (m > j) > 
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and e(2,fc) = 1, k = 1; = —1, k = 2. In (|7.8|) we have used the components of 
the determinant A(n, 2), 

(7.10) Z(l,y) = (l + 3y + 2y 2 ): Z(2, y) = (1 - y), 
and have written 

2 2 

(7.11) Z(fc,y) = ^^^,% i . 

i=0 k=l 

Consequently we may write 

2 4 np(m) 

(7.12) Z>(n,2) = ^^ £ C(*;m,j)l/ i /(m > j) 1 

i=0 m=0 j = l 

where 

2 

(7.13) C(i; m, j) = £ e(2, *;)*(», fc)T(m, j, k). 

fc=i 

To study the behavior of D(n, 2) for large n that is implied by f|T. 12[) . we note 
that the coefficients {C(i;m,j)} are independent of n, y = l/(2n), and I(m,j) 
is given by (|7.9|l . This integral is the product of two integrals each over a single 
variable of the form 

P oo 

(7.14) / d<$(t)t 2 "- 2 (t-r) r ', 

Jo 

with r\ being a small integer. For large t (|2.ip shows that <&(<) is very small, 
while t 2n is large, with the result that the integrand has a sharp maximum at 
t = t where 27rexp(4r). Using the Laplace method (see [I], [5]) we find that the 
largest terms in (|7.12j) occur when (i,m) = (1,0), (0,2). The ratio of the terms 
(l,0)/(0, 2) = 0(r),n — > oo. Calculations confirm this prediction, including the 
magnitude of the coefficients of these terms. 

It is important to note that C(0; 0, 1) = 0, as follows immediately from (IT. 13[) . If 
this coefficient were not zero, then the ratio (0, 0)/(l, 0) = 0(n), n —> oo. To obtain 
the correct behavior at very large n from (|7.13p it is essential that C(0; 0, 1) = 0. 

7.2. The method of Sec. 7.1 for order 2 can be extended to any r > 2 in a straight- 
forward manner. We repeat the main steps for the general case. We write 

(7.15) D(n, r) = y r{r -^ m g(i, n) r_i J T(2n + l) r D{n, r) 

Let q(i, k)), i = 1, 2, . . . , r be the permutation of tbe numbers 1, 2, . . . , r corre- 
sponding to index k = 1, 2, . . . , r!. The order in which these permutations are listed 
is immaterial. Define 

(7.16) u(i,k) = q(i,k) -i, i = l,2,...,r. 

We also define e(r, fc), k = 1, 2, . . . , r! to be equal to ±1 depending on whether the 
permutation is even (+1) or odd (—1). In place of (17.51) we choose 



r 

(7.17) R(k,t) = S]J ^ 2 M^)+r-i) 



fc = l,2,. 



i=i 
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where t = (t\, tz, . ■ ■ , t r ), and S means symmetrize. Then the components of 
det [bj~i+ n } i j=1 ^ r are 

(7.18) B(k) = S^p_jytMU)tf l - 2r+2 ^ R(k,t), fc = l,2,...,H 
For general r the analog of (|7.7I) is 

2r(r — 1) np(m) 

(7.19) R(k,t) = r 2r ^ Y,T{m,j,k)v(m,j), fc = 1,2, . . . ,H 

m= j— 1 

which defines the set of coefficients {T(m,j, k)} in this case. The determinant 
A(n,r) is a polynomial in y of degree r(r — 1), as stated in Lemma [6761 above, so 
for its components we have 

r(r-l) r< 

(7.20) Z(k,y)= J2 z ^ k ^ 1 - 

i=0 k=l 

Combining the above information just as in Sec. 7.1, we find that 

r(r — 1) 2r(r — 1) np(m) 

(7.21) D(n,r)= ]T £ £ C(i; m,j) y%m,j), 

i—0 m— j—1 

where 

r! 

(7.22) C(i; to, j) = ^ e(r, *;)*(*, k)T(m,j, k), 

fe=i 

and 

(7.23) I(m,j) = r 2r ^ ijlJ^dtMU)^- 2 ^ 

Calculations for small r suggest that the behavior noted for r = 2 may be 
generalized to all r. We propose 

CONJECTURE 3. For j = 1, . . . , np(m) the coefficients 

(7.24) C(i;m,j) = 0, m = 0, 1, . . . ,r(r - 1) - 2i - 2; i = 0, 1, . . . , r(r - l)/2 - 1; 

Remark 7.1. The specific values of the cofficients {& n } can have no influence on 
the proposed relation. Only the structure bj^i +n and the T function in (16.241 ) can 
matter. 

We predict that, given the truth of Conjecture 3, the Laplace approach men- 
tioned in Sec. 7.1 will lead to a rigorous demonstration that, for fixed order r, the 
determinant D(n,r) > for sufficiently high n. 

7.3. We expect that the upper representation will lead to a proof of Conjecture [2] 
This has already been demonstrate^ for some of the simpler cases. 



3 Will be available at http://publish.uwo.ca/~jnuttall/upper.pdf 
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Now we return briefly to the approach mentioned at the beginning of Sec. 6. In 
(|2.20|) we gave the relation between Ar p i(s,t) and the compound kernel derived from 
K(u, v) — $(tt + v), i.e. 



(7.25) A[p](s,t)=/ du dv(j)[ p }{u,^K [p] (u,v)(j3[ p }{v,£), 

J o Jo 

where u, v, s, t are defined as in (|2.18l) . In (|3.4p and (|3.8[> we gave examples showing 
how A[ p ](s,i) can be related to D(n,r) for r — 3,4. So long as n is large enough 
this procedure may be generalized to any value of order r. However we know that 
K(x, y) is not sign-regular for r > 4, so that (|7.25[) cannot be used directly to show 
that D(n, r) > for such r. 

Nevertheless we believe that it may be possible to use (|7.25[) to prove D(n, r) > 
so long as n is large enough for a given value of r. Just as in our application of the 
Laplace method we conjecture that, for large n, the dominant contribution to the 
integral in (|7.25|) arises from locations where all the components of (u, v) are large. 
Thus, for large n, it should be possible to find a quantity \i such that 



(7.26) A[pj(s,t;/i) = / du dv (j>[ p ](u, s)K [p] (u,v)(j) [p] (v,t), 

is a good approximation to A[ p ](s, t). Moreover, if fj, is large enough only a small 
error will be introduced in (|7.26j) by replacing $(w + v) in K{u, v) by the first term 
in (|2.1[) . When [i is large enough that the above approximations are adequate, then 
we can apply the notions of sign regularity to the reduced region of integration. It 
follows that D(n, r) > if Conjecture 3 holds, so that e r W r (y) > for large enough 
n, where W r {y) is given by (|6 . 2|) . 

8. DISCUSSION 

8.1. The results of this report, both rigorous and conjectured, are a testament to 
the power of computer-assisted mathematics. In several instances we find that, if 
we look beneath the surface, there are reasons why the behavior required by the 
RH occurs. The RH is perhaps not quite as mysterious as it often seems, but no 
doubt many secrets remain. 

The report raises a number of questions/problems that could led to further 
progress. We intend to list the problems at |http: / / publish.uwo.ca/ ~jnuttall where 



will also be found the unpublished documents listed in the References. 
The author would appreciate comments and corrections. 
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